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^Drainage  characteristics  of  grooved  runways  were  investigated.  An 
analytical  model,  developed  from  unsteady  flow  theory,  was  solved  on  a 
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30-ft-by-15-ft  Portland  cement  concrete  slab  sloping  at  1.5%  along  the  30-ft 
side.  Rainfall  rates  of  1.15  and  2.45  In. /hr  were  applied  to  the  slab  with 
rectangular  grooves  parallel  to  the  30-ft  side  and  spaced  at  •,  5,  2.5,  and 
1.25  In.  Data  on  the  slab  surface  showed  that  water  depths  were  smaller,  and 
the  hydraulic  roughness  was  larger,  than  expected.  The  model  was  adjusted 
accordingly,  and  used  to  generate  curves  of  water  depths  and  their  reductions 
versus  distance  out  to  100  ft  for  a  rainfall  of  3  In. /hr  and  a  hydraulic 
roughness  of  0.056.  Predicted  reductions  In  water  depths  ranged  from  100% 
through  39  ft  to  27%  at  100  ft  for  the  1.25- in.  spacing;  from  100%  through  19 
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at  100  ft  for  the  5-1 n.  spacing. 
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1.  INTRODUCTION 

1.1  PURPOSE  OF  THE  STUDY 

It  is  generally  known  that  the  wet  braking  action  of  an  aircraft  tire  is 
considerably  better  on  a  grooved  runway  than  on  an  ungrooved  runway. 

Aircraft  hydroplaning  is  eliminated  or  delayed  to  higher  speeds  on  runways 
grooved  normal  to  the  flight  path.  This  Improved  braking  performance  is 

i  * 

i  believed  to  be  the  result  of  a  dual  process  of  water  removal  from  the  tire¬ 

runway  Interface.  First,  the  grooves  Influence  the  surface  water  drainage  by 
providing  smooth  channels  Into  which  and  through  which  water  can  flow  freely. 
Whether  decreases  In  groove  spacing  will  continuously  Improve  the  drainage  on 
the  runways  Is  not  known.  The  smaller  groove  spacing  will  create  more  flow 
area  and  thus  more  flow  per  unit  length,  but  a  larger  spacing  might  be 
equally  effective  and  certainly  less  expensive.  Groove  spacing,  surface 
texture,  and  slope  of  the  runway  are  Interrelated  In  establishing  the  free 
flow  of  water.  Second,  the  grooves  provide  forced  water  escape  from  the 
tire-runway  interface  vrfien  the  aircraft  is  decelerated  on  a  water-covered 
runway.  Since  the  maximum  amount  of  water  that  can  be  removed  from  the 
runway  In  a  given  time  Is  limited,  both  the  free  flow  and  forced  escape  of 
water  are  Important. 


The  relative  Improvement  In  the  braking  performance  of  an  aircraft 
resulting  from  forced  water  escape  In  grooves  of  various  spaclngs  was 
established  In  Reference  1.  The  results  Indicated  that  an  Increase  In  groove 
spacing  causes  reduction  In  available  friction  level.  The  effect  of  groove 
spaclngs  on  the  free  flow  of  water  In  terms  of  the  water  depths  on  a  runway 
Is  the  subject  of  the  research  reported  here. 

The  research  was  divided  Into  two  phases.  In  Phase  I,  a  mathematical 
model  was  developed  which  Included  the  following  parameters:  (1)  transverse 
slope  of  surface,  (2)  surface  texture,  (3)  groove  size  and  shapes,  (4)  groove 
spacing,  and  (S)  uniform  rain  rate.  The  model  used  the  kinematic  wave 
approximation  to  the  shallow  water  equations,  for  grooved  and  ungrooved 
runways  with  a  zero  longitudinal  slope  and  a  transverse  slope  of  1.5X.  The 
assumption  was  made  that  all  rain  falling  on  the  pavement  drains  freely  from 
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the  edge  of  the  pavement,  except  for  water  that  either  is  stored  in  grooves 
or  merely  contributes  to  the  wetness  of  the  pavement.  Grooves  were  assumed 
to  run  through  the  edge  of  the  pavement;  however,  the  model  would  handle  the 
commonly  occurring  case  where  grooves  do  not  run  out  to  the  edge  of  the 
pavement.  Water  can  be  visualized  as  spilling  out  of  the  ends  of  these 
sloping  grooves  onto  the  ungrooved  portion  of  the  pavement.  Flow  can  then 
exit  from  the  pavement  edge  at  the  same  depth  it  would  have  on  an  ungrooved 
runway.  Numerical  solutions  of  the  model  were  generated  by  a  FORTRAN^  IV 
computer  program.  Hydroplaning  analysis  was  not  within  the  scope  of  this 
research,  even  though  the  results  reported  here  are  Important  to  that  topic. 

A  discussion  of  hydroplaning  may  be  found  In  publications  such  as  Agrawal  and 
Dalutolo  [1]  and  Horne  [2].*  Phase  I  has  been  documented  In  References  [3] 
and  [4], 

Phase  II  constituted  the  experimental  part  of  the  project.  This  phase 
Incorporated  the  following  elements: 

1.  Indoor  facility  using  a  rain-producing  machine  developed  In  a 
cooperating  study  sponsored  by  the  Federal  Highway  Administration, 
with  storm  durations  and  Intensities  being  restricted  to  those  In 
that  study. 

2.  Portland  cement  concrete  test  slab  approximately  15  ft 
(longitudinal)  by  30  ft  (transversal),  the  size  being  restricted  by 
the  rain  equipment  described  in  (1)  above. 

3.  Sampling  from  the  groove  spacing  and  rain  parameters  listed  for 
Phase  I,  for  a  fixed  transverse  slope  of  1.5X  and  an  as-built 
concrete  surface  with  0.25-1n. -square  grooves  only. 

4.  Determination  of  time  to  reach  steady  state  under  various  rain 
rates . 


•Numbers  In  brackets  refer  to  corresponding  items  In  the  list  of  references. 


5.  Measurement  of  Mater  depths  and  flow  rates  under  steady-state 
conditions. 

6.  Modification  of  the  Phase  I  mathematical  model  on  the  basis  of  the 
experimental  results  for  rectangular  grooves,  and  extrapolation 
beyond  the  range  of  experimental  parameters. 

Phase  II  has  been  documented  in  Reference  [5]. 

1.2  CONCEPTUAL  APPROACH  TO  THE  PHASE  I  ANALYTICAL  MODEL 

One-dimensional  flow  over  a  sloping,  ungrooved  section  of  pavement  is 
relatively  easy  to  visualize.  A  uniform  rainfall  begins.  Flow  and  depths 
over  the  whole  section  gradually  build  up  with  time  and  distance  downslope 
(unsteady  flow).  After  a  short  time,  steady-state  (equilibrium)  conditions 
are  reached  at  all  points  on  the  surface.  Application  of  the  equations  from 
the  kinematic  wave  approximation,  shown  later  In  detail,  occurs  smoothly, 
efficiently,  and  as  accurately  as  estimates  will  permit.  Maximum  values  of 
depth  and  flow  exist  at  the  downstream  end  of  the  pavement. 

Flow  over  the  grooved  pavement  Is  almost  as  easy  to  visualize,  with 
unsteady  conditions  followed  by  equilibrium  conditions.  Application  of  the 
one-dimensional  kinematic  wave  equations  does  not  occur  smoothly,  however, 
since  some  water  will  flow  laterally  Into  the  grooves.  This  lateral  inflow 
refers  to  local  flow  approximately  parallel  to  the  flight  path.  The  Inflow 
can  be  Imagined  to  occur  at  a  faster  rate  further  downstream  where  larger 
depths  (heads)  are  available,  and  will  add  to  the  flow  the  grooves  carry  due 
to  rain  falling  directly  into  them.  At  some  point  the  grooves  fill  up,  and 
flow  apparently  occurs  over  the  surface  In  a  manner  similar  to  that  on  the 
ungrooved  surface,  although  at  lower  depths. 

A  mathematical  model  for  the  grooved  pavement  was  attempted  by  using  a 
side-channel  weir  equation  to  estimate  the  lateral  Inflow.  Execution  of  the 
computer  program  for  this  model  appeared  successful  until  conservation  of 
mass  (continuity)  was  checked.  The  flow  leaving  the  pavement  was  less  than 
the  total  flow  occurring  on  the  pavement  due  to  rainfall.  After  several 
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refinements,  the  error  was  reduced  to  6%,  but  further  study  was  unsuccessful 
In  eliminating  the  error  and  the  model  was  set  aside. 

A  second  effort  was  attempted.  In  which  all  of  the  rain  falling  on  the 
pavement  at  the  upstream  end  was  compared  with  the  total  volume  of  the 
grooves  to  estimate  where  the  grooves  were  filled.  Once  this  point  was 
established,  computations  were  executed  beginning  there  as  though  the 
pavement  were  ungrooved,  except  that  the  flow  In  the  full  grooves  was 
accounted  for.  A  small  continuity  error  of  about  0.5X  was  found,  perhaps  due 
to  round-off  error.  This  model  was  set  aside  because  It  seemed  unrealistic. 
It  Is  difficult  to  Imagine  completely  filled  grooves  at  the  upstream  end  of 
the  pavement  with  no  flow  on  the  adjacent  textured  pavement.  Also,  at  the 
extreme  upstream  edge,  there  would  not  have  been  sufficient  rainwater  to  fill 
the  grooves. 

A  third  model  was  developed  using  the  concept  of  the  wetted  perimeter, 
l.e.,  the  total  linear  length  of  the  solid  boundary  of  a  cross  section 
through  which  flow  can  occur.  The  wetted  perimeter  Is  the  primary  factor 
related  to  the  boundary  shear  stress  caused  by  the  flow  of  water  over  the 
surface.  The  model  assumes  a  planar  ungrooved  boundary  which  has  a  width 
equivalent  to  the  wetted  perimeter  of  a  grooved  boundary,  thereby  preserving 
the  shear  area.  For  example.  In  Figure  1  (which  shows  the  sizes  and 
spaclngs,  s,  of  the  two  groove  patterns  considered  In  this  study,  as  well  as 
the  lateral  Inflow  Into  grooves  and  the  surface  depth),  the  width  of  the 
equivalent  planar  section  corresponding  to  s  for  the  rectangular  grooves  Is 
(s  +  0.5  In).  Thus,  a  segment  0.75  In.  In  width  models  the  groove  on  this 
planar  section.  Rain  still  falls  Into  the  middle  0.25  In.  of  this  segment, 
but  the  rain  falling  on  the  additional  0.25  In.  on  each  side  of  the  segment 
Is  also  considered  to  be  carried  In  the  groove  as  an  arbitrary  allotment  of 
lateral  Inflow  from  adjacent  surfaces.  Hence,  there  Is  available  water  to  be 
carried  on  the  adjacent  textured  surface,  but  less  than  on  an  ungrooved 
pavement.  A  computer  program  executed  this  model  successfully,  with 
continuity  exactly  satisfied. 

This  planar  model  Is  not  without  weaknesses.  The  allotment  of  lateral 
Inflow  Into  a  groove  from  adjacent  surfaces  Is  based  solely  on  how  much 


larger  is  the  wetted  perimeter  of  a  groove  compared  with  its  top  width.  The 
rectangular  groove  shown  in  Figure  1(a)  has  a  wetted  perimeter  three  times 
Its  top  width;  however,  the  wetted  perimeter  of  the  reflex-percussive  groove 
shown  in  Figure  1(b)  is  only  about  2556  larger  than  Its  top  width.  If  the 
region  of  Influence  of  a  groove  for  lateral  inflow  is  defined  as  the  wetted 
perimeter  minus  the  top  width,  then  the  region  of  influence  for  the 
rectangular  grooves  Is  0.5  In.,  whereas  for  the  reflex-percussive  groove  it 
Is  only  0.184  In.  One  would  think  that  the  two  should  be  equal,  but  the 
arbitrary  allotments  of  lateral  inflow  to  the  grooves'  shapes  will  be 
different.  Experimental  correction  of  the  mathematical  model,  however,  was 
rather  simple. 

There  are  two  other  minor  problems  with  the  planar  model  when  the 
model's  results  are  reverted  to  the  actual  situation  of  flow  over  the  grooved 
pavement.  The  first  Is  that  while  the  grooves  are  filling,  not  all  of  the 
wetted  perimeter  of  a  groove  Is  shear  area.  Second,  lateral  Inflow  to  a 
groove  will  continue  realistically  after  the  groove  Is  filled  up  and  until 
the  planar  model  depth  Is  equal  to  the  depth  above  a  groove;  but  beyond  that 
point,  lateral  Inflow  to  the  grooves  must  be  stopped. 

All  three  conceptual  models  assumed  a  constant  hydraulic  roughness 
coefficient  (Manning's  n)  across  a  cross  section,  even  though  the  rectangular 
grooves  will  be  shown  later  to  have  a  lesser  value  than  that  of  the  textured 
surface  as  a  result  of  polishing  while  being  grooved.  The  larger  Manning's  n 
value  was  considered  to  be  conservative  for  Phase  I  because  It  predicted 
larger  depths. 

The  model  based  on  the  concept  of  wetted  perimeter  was  chosen  for  both 
phases  of  this  research.  The  reasons  for  selecting  it  are  the  ease  with 
which  continuity  Is  satisfied  and  the  facility  with  which  adjustments  can  be 
made  to  It  based  on  experimental  evidence. 

1.3  SCOPE  OF  THE  PHASE  II  EXPERIMENTS 

The  project  called  for  the  construction  of  an  Indoor  Portland  cement 
concrete  slab  and  the  provision  of  all  the  necessary  equipment  for  runoff 


measurements.  The  slab  was  built  at  the  Pennsylvania  Transportation  Research 
Facilities.  Rainfall  simulation  equipment,  developed  in  a  cooperating  study 
for  an  area  30  ft  by  15  ft  to  provide  a  rainfall  Intensity  of  about  1  inch 
per  hour,  supplied  the  artificial  rain.  The  slab  size  of  30  ft  by  15  ft  was 
a  cost-effective  accommodation  to  the  rain  simulation  equipment,  since  the 
30-ft  dimension  Is  not  typical  of  the  100-ft  transverse  width  of  a  runway 
lane. 


The  complete  project  work  plan  consisted  of  two  rainfall  rates 
(approximately  1  In. /hr  and  2.5  In. /hr)  applied  to  each  of  the  four 
rectangular  groove  spaclngs  shown  In  Figure  1(a).  Water  depths  were  recorded 
at  predetermined  locations  on  the  slab  where  sensor  openings  were  connected 
to  pressure  transducers.  FAA  personnel  cut  the  grooves,  with  their 
equipment,  parallel  to  the  30-ft  dimension.  No  experiments  were  to  be 
conducted  on  the  ref lex -percussive  grooves  shown  in  Figure  1(b),  although 
prediction  of  their  effect  on  runoff  was  to  be  made  from  the  experience  on 
rectangular  grooves.  The  spacing  of  1.25  In.  has  been  widely  used  and  formed 
the  basis  for  the  spaclngs  studied.  See  Reference  [6]  for  a  discussion  of  3- 
1n.  spacing  versus  1.25- in.  spacing. 

1.4  RELATION  OF  PHASE  II  TO  PHASE  I 


It  was  anticipated  that  the  experiments  would  yield  Information 
concerning  the  amount  of  lateral  Inflow  to  grooves,  and  the  hydraulic 
roughness  coefficient  for  pavement  and  groove  surfaces.  A  computer  program 
developed  to  solve  the  analytical  model  was  modified  to  accommodate  this 
Information.  In  addition,  the  modified  model  was  used  to  generate 
mathematical  water  surface  profiles  to  be  compared  with  plotted  data  points. 
Finally,  the  modified  model  was  used  to  predict  runoff  depths  beyond  the  30- 
ft  range  of  the  slab  for  a  variety  of  conditions  of  pavement  roughness  and 
rainfall  rate  for  the  two  groove  shapes  shown  In  Figure  1. 


2.  ANALYTICAL  MOOEL 


2.1  BACKGROUND 

The  equations  for  predicting  water  thickness  on  an  overland  flow  surface 
such  as  a  sloping  pavement  are  known  as  the  shallow  water  equations.  These 
equations  are  based  on  conservation  of  mass  and  momentun  during  unsteady, 
spatially  varied  flow,  and  have  the  one-dimensional  form  shown  In  Henderson 
[7]  and  Chow  [8]: 

Continuity:  ft  +  ulj+y-Jj-q  (1) 

Momentum:  |f+ul7+9|x“  9^So  *  V  “  ^  (u  "  v)  (2) 

where 

u  *  local  velocity  In  fps 
y  *  local  depth  In  ft 
g  •  gravitational  acceleration  constant 
q  ■  lateral  Inflow  or  rainfall  excess  rate  in  cfs/ft2 
S0  ■  slope  of  plane  or  channel  in  ft/ft 
Sf  ■  friction  slope  In  ft/ft 
v  ■  velocity  of  lateral  Inflow  In  ft/s 
x  *  distance  downstream  In  ft 
t  ■  time  In  s 

A  definition  sketch  of  the  variables  Is  shown  In  Figure  2.  Note  that  x  and  y 
are  referenced  to  a  horizontal  line  and  not  to  the  sloping  line  of  the  plane 
surface.  The  difference  Is  negligible  for  small  values  of  S0.  Equations  (1) 
and  (2)  have  been  applied  to  many  problems  of  open-channel  flow.  Including 
the  analysis  of  flood  wave  movement  In  river  systems.  Their  solution  on  the 
computer  requires  the  use  of  complex  finite-differencing  methods  which  are 
time-consuming  and  sometimes  numerically  unstable,  as  Liggett  and  Woolhlser 
[9]  have  Indicated. 


As  will  be  shown  in  Section  2.2,  many  Investigators  have  examined  the 
properties  of  the  kinematic  wave  solution  as  an  approximation  to  the  complete 
shallow  water  equations.  Typical  results  are  given  by  Llghthi 1 1  and  Whitham 

[10] .  The  kinematic  wave  approximation  is  generally  valid  when  the  friction 
forces  on  the  plane  or  channel  are  assumed  to  be  just  balanced  by  the 
gravitational  forces.  Hence,  S0  ■  $f,  and  vrfien  the  slopes  disappear  from 
the  momentun  equation,  the  resulting  absence  of  force  terms  Implies 
"kinematic."  The  ever-changing  y  In  Figure  2,  caused  by  rainfall  excess. 
Implies  "wave."  Under  these  conditions,  the  Chezy  or  Manning  formula 
(described  later  as  equations  (6)  and  (7))  can  be  used  to  describe  the 
resulting  flow  if  it  Is  fully  rough  turbulent,  as  justified  by  the  condition 
In  equation  (11).  This  situation  Is  usually  the  case.  Woolhlser  and  Liggett 

[11]  In  particular  have  tested  kinematic  solutions  for  the  rising  hydrograph 
on  an  overland  flow  plane  and  have  concluded  that  most  cases  of  overland  flow 
can  be  represented  by  the  mathematically  simpler  kinematic  wave 
approximation. 

Numerous  studies  of  the  kinematic  wave  equation  have  been  completed 
since  the  late  1960s,  and  this  simplified  approach  has  been  almost 
universally  accepted  for  the  analysis  of  surface  runoff  and  overland  flow 
from  paved  surfaces.  Klbler  and  Aron  [12]  have  applied  a  kinematic  wave 
runoff  model  to  the  controlled  laboratory  runoff  studies  of  Izzard  [13]  for 
various  pavement  surfaces.  It  can  be  seen  from  the  typical  results  of  the 
kinematic  wave  model  shown  In  Figure  3  that  simulated  runoff  agrees  closely 
with  the  measured  Izzard  runoff  under  changing  rainfall  supply  rates  at  two 
different  values  of  hydraulic  roughness  (Manning's  n). 

Gallaway,  Schiller,  and  Rose  [14]  have  correctly  Identified  the  pavement 
runoff  factors  as  rainfall  Intensity,  pavement  texture,  cross  slope,  and  flow 
length,  and  have  Incorporated  these  into  an  empirical  expression  for  average 
water  depth  above  the  top  of  the  texture.  Although  their  study  Is  valuable 
as  a  background  for  the  present  study.  It  does  not  account  for  pavements  that 
are  grooved,  nor  does  It  Indicate  a  basic  measure  of  flow  resistance,  such  as 
Manning's  n,  which  Is  vital  to  a  kinematic  wave  solution.  Reed  and  Klbler 
[15]  have  used  data  from  this  and  other  studies  [16,  17]  to  calculate 
Manning's  n  using  the  equations  shown  In  Section  2.2,  and  subsequently  to 


Figure  3.  Kinematic  Wave  Runoff  Model  Application  to  Izzard  Asphalt  Surface  [12) 
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Investigate  Its  prediction.  Some  of  their  results  are  shown  In  Figure  4. 

The  average  and  upper  limit  curves  of  Manning's  n  values  were  used  In  this 
study  to  estimate  hydraulic  roughness.  The  texture  depth  can  be  obtained 
from  macrotexture  measurements,  which  are  typically  described  In  Rose  et  al. 
[18]  and  Henry  and  Hegmon  [19].  Representative  values  of  Manning's  n  were 
selected  Initially  for  the  analytical  model  from  typical  runway  macrotextures 
between  0.01  In.  and  0.03  In;  however,  measurement  of  the  actual 
macrotexture  of  the  experimental  slab  In  Phase  II  was  necessary  to  obtain  n 
from  Figure  4  In  the  early  experimental  analysis. 

In  their  study,  Reed  and  Klbler  [15]  have  used  conditions  described  In 
the  next  section  by  equations  (3)  and  (11)  to  justify  their  use  of  the 
kinematic  wave  approximation  as  well  as  the  use  of  the  Manning  equation. 

Even  In  cases  where  equation  (11)  might  not  be  satisfied,  the  flow  might 
still  be  turbulent  because  of  the  turbulence  added  by  pelting  rain,  a  view 
supported  by  the  Army  Corps  of  Engineers  [20}. 

2.2  BASIC  EQUATIONS  FOR  OVERLAND  FLOW  ON  PAVED  SURFACES 

When  rain  falls  on  pavements  such  as  highways,  streets,  and  runways, 
flow  over  the  surfaces  occurs  In  a  form  that  has  become  known  as  urban 
runoff.  This  flow  Is  analogous  to  the  flow  that  would  occur  In  an  Impervious 
open  channel  of  very  great  width,  and  hence  may  be  referred  to  as  overland 
flow.  The  hydraulics  of  overland  flow  are  generally  described  by  the 
kinematic  approximation  to  equations  (1)  and  (2).  It  has  been  shown  [11] 
that  the  kinematic  wave  approximation  to  these  equations  Is  generally  valid 
when  the  parameter  K  exceeds  20  and  the  Froude  number,  F0,  exceeds  0.5, 
where 

K  .  so.--L-a  (3) 

H0  *  Fq2 

where 

$0  ■  slope  of  flow  path,  ft/ft 
l0  ■  length,  ft 
H0  >  normal  depth,  ft 
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F0  ■  Froude  number  «  V0/  /gH0 
V0  ■  normal  flow  velocity,  ft/s  ■  (q  •  L0)/H0 
q  ■  lateral  Inflow  or  rainfall  excess  rate  In  ft/s 
g  ■  gravitational  constant  ■  32.2  ft/s2 

Under  subcrltlcal  flow  conditions  (F0  <  1)  where  F0  Is  less  than  0.5,  an 
additional  criterion  must  be  met  [21]: 


F02  K  •  S0  L0/H0  >  5 


Except  on  extremely  flat  slopes  with  Intense  rainfall,  most  cases  of 
overland  flow  on  paved  surfaces,  Including  the  ones  In  this  study,  will  fall 
easily  Into  the  kinematic  range  based  on  the  preceding  criteria.  This  means 
that  the  more  complex  equations  for  spatially  varied,  unsteady  flow  can  be 
replaced  by  the  simple  equations  Involving  mass  continuity  and  uniform  flow. 
The  partial  differential  form  of  the  kinematic  equation  Is  well  known  and  has 
been  used  frequently  for  overland  flow  modeling  In  urban  watersheds  [22]. 
Integrating  the  characteristic  equation  In  an  x-t  plane  under  constant 
rainfall  condition  yields  a  simple  relation  for  the  equilibria  flow  depth 
(Yeq)  at  different  points  In  the  downstream  direction: 


Y  *  q  t 
eq  M  eq 


where 


q  ■  lateral  Inflow  or  rainfall  excess  rate  In  cfs/ft2 
teq  *  equilibrium  flow  time,  s 

The  equilibrium  flow  time,  t^q,  on  an  Impervious  surface  under  constant 
rainfall  excess  rate,  q.  Is  determined  by  the  equation: 


where 


L  ■  length  along  flow  path  or  width  of  runway,  ft 
a,  m  ■  parameters  of  the  hydraulic  friction  relationship 
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Under  kinematic  flow  conditions,  the  friction  or  energy  slope,  Sf»  is 
equal  to  the  slope  of  the  plane,  S0,  and  this  permits  the  use  of  a  simplified 
uniform  flow  velocity  relation  of  the  form: 


V  - 


(6) 


where 

V  «  normal  flow  velocity,  ft/s 

Y  ■  flow  depth,  ft 

For  uniform  turbulent  flow  on  an  overland  surface,  the  parameters  o  and  m  can 
be  determined  from  the  Manning  equation  as 

a  .  Lii  s0°*5  (English  units)  (7) 

n 

m  *  5/3  (dimensionless  constant) 

Equations  (4),  (5),  (6),  and  (7)  provide  the  theoretical  basis  for  the  flow 
model  described  earlier  In  Section  1.2. 


Under  the  assumption  of  turbulent  flow,  equation  (7)  can  be  applied  to 
equation  (5)  for  teq  to  obtain 


teq(sec) 


56.25  L 


1 


0.4 


0.6n0.6 

c  0.3 
^0 


(8) 


where 

1  ■  rainfall  Intensity  In  In. /hr 
n  *  Manning  roughness  coefficient 

and  the  other  terms  are  as  defined  earlier.  Equation  (4)  for  Yeq  can  now 
be  written 
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Discharges  under  uniform  rain  can  now  be  computed  for  all  points  L  ft 
downstream  using 


Q  .  o  y"  (l0) 

where  Q  •  discharge,  cfs/ft,  and  the  other  terms  are  as  defined  earlier. 
Equilibrium  flows  can  also  be  obtained  from  continuity  considerations  by 
calculating  the  Inflow  onto  the  pavement  upstream  from  Intensity.  Equations 
(8),  (9),  and  (10),  together  with  (6)  and  (7),  provide  the  computational 
basis  for  the  FORTRAN  IV  program  written  to  execute  the  flow  model. 

Henderson  [7]  has  suggested  that  fully  rough  flow  conditions  prevail 

when 

n6  /RhSf  >  1.9  x  10~13*  {11) 


where 

Rh  ■  hydraulic  radius  (assumed  equal  to  overland  flow  depth,  Y) 

Sf  *  friction  slope  (Sf  ■  SQ)  under  kinematic  flow  conditions 

Equation  (11)  Is  satisfied  In  the  present  study  for  most  of  the  data, 
allowing  for  some  reasonable  minimum  depth  such  as  0.02  In.  Hence,  the  use 
of  the  Manning  equation  to  describe  the  friction  relations  for  the  pavements 
In  this  study  Is  justified.  For  Henderson  numbers  less  than  1.9  x  10~13, 
the  flow  Is  transitional  or  laminar,  and  consequently  Manning's  n  should  not 
be  used  unless  It  Includes  the  Reynolds  number  [23]. 

2.3  COMPUTER  SOLUTIONS  OF  THE  ORIGINAL  MODEL 


A  FORTRAN  IV  program  was  written  to  solve  the  mathematical  flow  model. 
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A. 


The  program  was  executed  on  an  IBM  370/Model  3033  system  to  generate  flow 
profiles  which  can  be  depicted  graphically.  Several  graphs  were  produced  for 
a  wide  range  of  values  of  the  pertinent  parameters.  These  graphs  are  shown 
In  Reference  [3],  but  they  are  not  completely  reproduced  here  because  of  the 
eventual  modification  of  the  model  from  experimental  experience,  particularly 
with  grooved  pavement.  However,  two  of  the  graphs  are  shown  in  Figures  5  and 
6.  These  curves  show  that  Increasing  rainfall  rate  and  Manning's  n  values 
Increase  water  depths  on  ungrooved  runways.  The  curves  are  valid  for  the 
values  of  the  parameters  shown,  despite  later  modifications  to  the  model. 

The  figures  were  plotted  using  a  Hewlett-Packard  graphics  terminal  which 
selects  appropriate  scales  for  the  axes.  The  scales  are  not  always 
convenient  from  the  standpoint  of  good  plotting  practice.  Axes  are  labeled 
In  word  descriptions  for  clarity. 

Note  that  the  maximum  rainfall  rate  of  4  in. /hr  In  Figure  5  was  not  the 
maximum  considered  In  plotting  some  of  the  other  graphs  (not  shown).  That 
maximum  was  6  In./hr,  which  was  chosen  because,  together  with  the  length  of 
the  flow  path  (which  Is  related  to  rainfall  duration),  it  statistically 
represents  a  storm  with  a  return  period  of  about  25  years  In  Pennsylvania. 
Such  a  storm  Is  typically  used  to  estimate  design  flows  for  storm  drainage 
systems.  Unfortunately,  the  rain  equipment  used  in  the  experiments  In  Phase 
II  has  a  limiting  output  of  about  3  In. /hr.  Hence,  that  value  was  used  for 
predictive  analysis  later. 


Effect  of  Variable  Rainfall  Rate  on  Water  Depth,  Y  ,  for  Ungrooved  Surface 


Effect  of  Variable  Manning  Roughness,  n,  on  Water  Depth,  Y  ,  for  Ungrooved  Surface 


3.  EXPERIMENTAL  SYSTEM 


3.1  CONCRETE  SLAB  AND  ARTIFICIAL  RAINFALL 


The  concrete  slab  was  constructed  to  be  4  In.  thick  with  light 
reinforcing  to  minimize  tension  cracks.  A  plan  view  Is  shown  In  Figure  7. 

The  slab  slopes  along  the  30-ft  dimension  at  1.5%.  The  slope  was  set 
accurately  on  the  concrete  forms  prior  to  construction,  using  precise 
surveying  equipment.  It  was  checked  again  after  the  concrete  had  been 
poured,  typically  finished  with  a  stiff  broom  in  the  30-ft  direction,  and  had 
hardened.  The  slab  was  set  on  a  raised  pad  of  crushed  limestone  aggregate 
embedded  with  polyvinyl  chorlde  (PVC)  pipes  as  shown,  with  the  top  of  pipes  at 
the  aggregate  pad  surface.  The  black  dots  shown  over  the  PVC  pipes  In  Figure 
7  are  3/8-in.  holes  drilled  through  the  slab  to  the  empty  PVC  pipes  to 
accommodate  Instrumentation  tubing  beneath  the  slab.  The  weir  box  was  used  to 
measure  the  flow  off  the  downstream  edge  of  the  pavement.  When  this  flow  was 
compared  with  the  Inflow  through  the  rainfall  simulator,  an  estimate  of  the 
loss  of  runoff  due  to  overspray  or  flow  off  the  side  edges  of  the  slab  could 
be  made.  This  comparison  is  not  reported  here,  since  it  is  not  considered 
relevant  to  this  study. 

Figure  8(a)  shows  a  sand-patch  test  being  conducted  on  the  finished 
ungrooved  surface  of  the  slab  to  determine  an  average  texture  depth  at  a 
given  location.  Fifteen  such  tests  were  performed  uniformly  on  the  surface 
of  the  slab,  and  the  results  showed  an  average  texture  depth  of  0.026  in., 
with  a  range  of  0.022  In.  to  0.029  In.  The  texture  depth.  In  conjunction  with 
Figure  4,  provides  an  estimate  of  the  hydraulic  roughness  coefficient.  A 
special  sand-patch  test  was  made  near  the  center  of  the  slab,  where  a  rather 
severe  blemish  occurred  In  the  broom  finishing  of  the  fresh  concrete.  The 
test  showed  the  blemish  to  have  an  average  texture  depth  of  0.048  in.  The 
blemish  runs  throughout  the  15-ft  dimension  of  the  slab  and  can  be  observed  in 
the  photograph  In  Figure  13(a),  showi  later  In  this  report. 

The  photograph  In  Figure  8(b)  shows  an  overall  view  of  the  experimental 
system  in  operation,  with  artificial  rain  falling  on  the  grooved  slab.  Eight 
tripods  support  the  spray  nozzles,  valves,  and  pressure  gages  of  the  rainfall 
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Figures.  Photographs  of  the  Concrete  Test  Slab  and 
Artificial  Rain  Equipment 
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simulation  equipment.  A  600-gallon  tank  truck  supplies  Mater  to  the  system. 
Tank  capacity  can  be  a  limiting  factor  In  the  duration  of  testing.  At  the 
extreme  left  center  of  the  photograph*  a  tripod  may  be  observed  standing 
partially  on  the  slab.  This  condition  was  created  by  an  error  In  the  slab 
location.  The  slab  was  not  placed  far  enough  away  from  the  wall  to  provide 
Space  for  setting  up  a  tripod.  The  error  posed  no  difficulty  In  the 
experiments  except  for  voiding  the  use  of  the  column  of  sensor  openings  marked 
NX"  In  Figure  7.  Also*  fewer  grooves  were  cut  on  that  side  of  the  slab. 

Figure  9  shows  the  grooves  being  cut  by  FAA  personnel  using  a  diamond -tipped 
saw  machine. 

The  rainfall  simulation  equipment  has  been  documented  by  Davis  [24].  To 
establish  an  Intensity  of  1  In. /hr,  four  nozzles  in  a  predetermined  pattern 
were  used,  and  the  pressure  gages  at  the  nozzles  were  regulated  to  12  psl. 

For  2.5  In. /hr,  eight  nozzles  were  used,  and  the  pressure  gages  were  set  at  20 
psl.  It  should  be  understood  that  these  rainfall  rates  are  spatial  averages 
and  not  absolutely  uniform.  Standard  deviations  In  these  averages  could  be 
about  ±0.3  In. /hr.  To  ensure  a  closer  estimate  of  the  rainfall  Intensities  in 
this  research,  collection  cup  data  similar  to  those  taken  by  Davis  were  taken 
in  the  Important  flow  regions  spanned  by  the  sensor  holes  on  the  slab.  These 
cups  may  be  seen  In  Figure  10(a),  along  with  a  sensor  opening  shown  at  the  tip 
of  the  pen  on  the  textured  slab  surface.  The  collection  cup  data  more 
accurately  defined  rain  rates  for  this  study.  These  rates  were  1.15  In. /hr 
and  2.45  In. /hr.  Standard  deviations  were  not  determined  for  these  averages, 
but  experience  has  shown  that  cup  data  collected  In  the  center  of  a  ralned-on 
area  are  more  nearly  uniform.  Figure  10(b)  Is  a  photograph  of  the  broom- 
finished  surface,  showing  a  rectangular  groove  and  a  more  discernible  3/8-1n. 
sensor  opening. 

3.2  PRESSURE  TRANSDUCERS 


Figure  11(a)  shows  the  flexible  tubing  connected  to  the  sensor  openings, 
emerging  from  beneath  the  slab  In  a  4- In. -diameter  PVC  pipe.  Figure  11(b) 
shows  one  of  the  tubes  being  connected  to  a  pressure  transducer.  Ten  of 
these  U.S.-made  transducers  were  purchased  for  this  study  because  their  cost 
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(a)  Collection  cups  and  sensor  opening 


(b)  Rectangular  groove  and  sensor  opening 


Figure  10.  Photographs  of  Textured  Surface  of  Slab  Showing 
Collection  Cups,  Sensor  Openings,  and  Rectangular 
Groove 
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was  comparable,  without  sacrifice  of  accuracy,  to  the  single  one  made  In 
Europe  and  proposed  Initially  for  these  experiments.  In  Figure  8(b),  the 
wooden  box  In  which  the  transducers  are  mounted  can  be  seen  (left  center) 
resting  on  the  floor  at  the  side  of  the  slab.  All  the  transducers  were  placed 
below  the  slab  surface  to  register  positive  pressures,  since  they  were 
calibrated  for  positive  pressures  In  this  study.  The  transducers  sense  small 
changes  of  pressure  In  the  tubing  (flow  depths  on  the  surface  of  the  slab)  and 
send  an  electrical  signal  to  a  digital  multimeter.  A  5  volt  X  power  supply 
device  in  the  form  of  a  function  generator  was  used  with  the  system.  The 
multimeter  was  a  Oata  Precision  Model  258,  with  an  LCO  readout  to  the  nearest 
0.01  mV. 

Appendix  A  contains  technical  data  on  the  pressure  transducers. 

Including  the  calibration  coefficients  (K)  determined  Individually  for  the 
transducers  In  this  study.  Figure  12  shows  the  linear  nature  of  the  data, 
which  determines  the  maximum  and  minimum  K  values  (slopes).  The  transducers 
were  calibrated  using  a  container  with  a  bottom  opening  to  which  a  transducer 
could  be  connected.  Water  levels  In  the  container  were  varied,  and  changes  in 
voltage  response  by  a  transducer  were  recorded  against  changes  In  water  depth 
measured  by  a  mechanical  point  gage  precise  to  0.001  ft. 

The  depth-measuring  system  using  pressure  transducers  Is  considered 
superior  to  any  system  that  1$  mounted  either  above  or  on  the  slab  surface. 
Surface  systems  can  obstruct  the  very  quantity  they  are  Intended  to  measure. 
One  of  the  disadvantages  of  transducers  Is  that  their  sensor  positions  are 
fixed.  Reference  [14]  describes  an  excellent  above-the-surface  system 
Involving  a  point  gage  with  a  variable  (x,y)  position. 

3.3  EXPERIMENTAL  PRXEOURE 

Two  types  of  experiments  were  conducted:  quantitative  experiments  which 
Involved  the  recording  of  numerical  data,  and  qualitative  ones  which  involved 
the  observation  of  water  movement  patterns  on  the  slab,  made  by  a  fluorescein 
dye.  Some  quantitative  work  also  involved  dye. 


(a)  Flexible  tubing  emerging  from  PVC  pipe 


(b)  Pressure  transducer  connected  to  sensor  tubing 

Figure  11.  Photographs  of  Connection  of  Pressure  Transducers 
to  Sensor  Openings  Using  Flexible  Tubing 
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3.3.1  Equilibrium  Water  Depth  and  Hydraulic  Roughness  Data 

Prior  to  activation  of  the  artificial  rainfall,  the  tubes  connecting  the 
transducers  to  the  slab  sensor  openings  were  checked  to  see  whether  they  were 
filled  with  water  and  had  no  air  pockets.  This  procedure  is  difficult  because 
much  of  the  tubing  is  concealed  in  the  PVC  pipes.  It  is  even  more  difficult 
for  the  upstrean  row  of  sensors  (7.5  ft  in  Figure  7),  where  depths  are 
extremely  small  and  as  much  as  25  ft  of  tubing  Is  connected  to  the  sensors. 

The  water  level  In  the'sensor  openings  in  the  slab  was  then  checked,  and  a 
judgment  was  made  as  to  when  the  level  was  at  the  average  texture  depth  datum. 
In  this  process  a  wetting  agent  In  the  form  of  a  mild  detergent  was  placed  on 
the  opening  with  an  eye  dropper  In  order  to  break  down  surface  tension  and 
reduce  capillarity  effects.  When  the  water  levels  were  judged  to  be  at  the 
proper  datum,  the  transducers  were  zeroed  for  this  level  by  recording  their 
output  voltage.  At  this  point,  the  rainfall  equipment  was  prepared  for  an 
experiment. 

With  the  upstream  valves  of  all  eight  nozzles  of  the  rainfall  simulation 
equipment  closed,  the  centrifugal  pump  on  the  tank  truck  was  started,  which 
applied  pressure  to  the  system  supplying  water  to  the  nozzles.  The  pipes  of 
translucent  material  were  checked  for  air  pockets,  which  were  then  purged 
through  small  bleed  valves  with  a  minimum  loss  of  water.  Then,  each  nozzle 
was  set  for  the  appropriate  pressure  by  opening  the  valve,  adjusting  the 
pressure  regulator,  and  closing  the  valve.  Again,  the  loss  of  water  is 
minimized  because  this  operation  takes  only  a  few  seconds.  When  all  the 
nozzles  were  set,  the  pump  was  turned  off.  The  valves  to  all  of  the  nozzles 
were  opened,  and  a  small  residual  spray  occurred  for  a  short  time  because  of 
the  release  of  pressure.  Testing  began  when  the  pump  was  started  and  the 
required  Intensity  sprayed  forth.  For  the  tests  reported  here,  the  truck 
capacity  provided  adequate  durations  of  rainfall,  particularly  with  the 
conservative  startup  procedure.  A  fully  running  system  could  purge  itself  of 

air  as  well  as  have  Its  pressures  set;  however,  the  amount  of  water  In  the 
tank  when  testing  began  might  be  very  inadequate. 

Equilibrium  water  depth  data  are  given  in  Appendix  B  for  the  ungrooved 
surface  and  the  grooved  surfaces.  (Voltage  drops  were  converted  to  depths 
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prior  to  being  recorded  in  the  tables;  see  Appendix  A.)  The  data  were 
recorded  as  average  depths  at  locations  indicated  in  Figure  7.  A  range  of 
data  values  is  shown  for  each  depth  and  was  determined  by  three  or  more 
Independent  depth  measurements  over  several  runs.  Data  were  not  recorded  for 
the  upstream  row  of  sensors  at  7.5  ft.  There  was  no  discernible  depth  for  the 
grooved  pavement  there,  and,  although  some  depth  could  be  sensed  for  the 
ungrooved  pavement,  the  readings  were  considered  unreliable.  Moreover,  there 
are  no  data  shown  In  Appendix  B  for  grooves  spaced  at  2.5  In.  under  rainfall 
of  1.15  In. /hr,  or  for  grooves  spaced  at  1.25  in.  under  rainfall  of  1.15 
in. /hr  or  2.45  In. /hr.  Repeated  measurements  for  these  conditions  resulted  In 
no  measurable  water  depths  anywhere  on  the  slab  surface. 

The  measured  water  depths  for  the  ungrooved  pavement  appear  appreciably 
larger  than  the  depths  obtained  from  the  mathematical  model  of  Phase  I.  A 
possible  reason  for  this  discrepancy  Is  the  uncertainty  in  the  value  of  the 
hydraulic  resistance  coefficient,  Manning's  n,  used  in  the  mathematical  model. 
If  the  average  texture  depth  of  0.026  in.  is  used.  Manning's  n  from  Figure  4 
is  0.042  for  the  median  line  and  0.056  for  the  upper  limit  line.  If  the 
texture  depth  of  0.048  In.  Is  used  for  the  finishing  blemish  mentioned 
earlier  (see  Figure  13a),  the  upper  limit  line  of  the  curve  yields  n  ■  0.081. 
Of  course,  the  blemish  occurs  only  In  one  area  of  the  pavement,  but  It  extends 
all  the  way  across  the  slab.  It  was  decided,  therefore,  to  make  a  dye  trace 
determination  of  Manning's  n  prior  to  grooving  the  pavement.  For  this 
determination,  the  time  to  reach  equilibrium  conditions  at  the  downstream  edge 
of  the  slab  Is  customarily  approximated  by  the  travel  time  of  the  dye  over  the 
entire  length  of  the  slab  under  equilibrium  conditions.  Manning's  n  was 
obtained  by  solving  equation  (8).  The  results  are  shown  In  Table  1.  The  n 
value  of  0.096  seems  high,  but  Reference  [25]  indicates  that  n  might  be 
between  0.05  and  0.15  for  water  depths  comparable  to  those  reported  here. 

Also,  In  a  separate  and  as  yet  unpublished  thesis  study  not  involving  dye  or 
rain,  n  varied  from  0.03  to  0.10  at  a  fixed  location  on  the  slab  surface. 
Hence,  the  n  values  used  here  are  compatible  with  those  reported  elsewhere. 

Not  suprlslngly,  the  water  depths  shown  in  Appendix  B  are  lower  for 
grooved  pavement  than  for  the  ungrooved  pavement.  However,  there  are  no 
surface  depths  upstream,  a  result  which  is  in  disagreement  with  the  unadjusted 
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Table  1.  Dye  Test  Results:  Determination  of  Manning's  n  for 

Ungrooved  Surface 


(Distance  =  30 

ft;  slope  * 

0.015) 

Rainfall  Intensity:  1.15  in. /hr 

Rainfall  Intensity:  2.45  in. /hr 

Trial 

Equilibrium  Time  (sec) 

Trial 

Equilibrltin  Time  (sec) 

1 

347 

1 

255 

2 

362 

2 

259 

3 

349 

Z 

267 

Average 

353 

Average 

260 

n  -  0.096 

n  *  0.096 

Table  2. 

Test  Results:  Determination  of  Manning's  n  for  Grooves 

(Distance  *  15  ft; 

grooves  flowing  full) 

Trial 

Travel  Time  (sec) 

1 

21 

2 

23 

3 

24 

4 

23 

5 

20 

6 

21 

Average 

22 

ng  *  0.0097  ■  0.010 
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mathematical  model.  This  condition  suggested  that  qualitative  dye  tests  be 
conducted  to  confirm  that  all  the  rain  upstream  was  being  carried  in  grooves, 
and  that  quantitative  dye  tests  be  run  to  determine  Manning's  n  for  grooves. 
The  qualitative  tests  are  discussed  In  the  next  section. 

A  small  study  was  undertaken  to  determine  n  values  for  groove  surfaces, 
which  are  polished  to  some  degree  by  the  cutting  process.  A  flow  was  created 
from  a  hose  laid  on  the  surface  of  the  grooved  slab  which  was  not  subject  to 
rainfall.  The  flow  was  adjusted  so  that  reasonably  long  stretches  of  the 
grooves  would  flow  full  without  overflowing  onto  the  surface.  The  Manning 
equation,  which  Is  represented  by  equations  (6)  and  (7)  for  a  modeled  groove 
In  the  Phase  I  model,  can  be  solved  for  n  If  the  flow  velocity  can  be 
measured.  Velocities  were  measured  by  timing  the  movement  of  a  tiny  piece  of 
styrofoam  moving  In  the  grooves  over  a  distance  of  15  ft.  The  results  of 
these  measurements  are  shown  In  Table  2.  The  ng  value  of  about  0.010  is  not 
surprising,  since  a  table  look-up  for  such  a  surface  might  show  a  value  of 
that  magnitude. 

3.3.2  Qualitative  Tests  Using  Dye 

A  series  of  dye  tests  was  undertaken  to  observe  water  movement  on  the 
grooved  slab.  The  dye  used  was  a  dark-red  powder  vrtilch  becomes  greenish- 
yellow  when  dissolved  in  water.  The  dye  may  be  used  in  dry  or  liquid  form, 
but  the  most  effective  observations  In  this  study  were  made  by  placing  the 
powder  along  the  upper  edge  of  the  slab  and  allowing  the  rainfall  to  pick  it 
up  and  diffuse  It  downstream. 


Two  of  the  color  photographs  taken  during  these  tests  are  shown  in  Figure 
13.  The  left  side  of  Figure  13{a)  shows  the  dye  being  dispersed  at  the 
upstream  end  of  the  slab.  Very  little  dye  is  visible  on  the  surface  in  the 
middle  of  the  photograph,  indicating  that  the  water  is  being  carried  in  the 
grooves.  On  the  right  side,  dye  appears  to  be  coming  out  of  the  grooves  onto 
the  surface,  vrfilch  Indicates  that  the  grooves  are  filled  and  overflowing.  The 
surface  blemish  Is  evident  In  Figure  13(a)  and  Is  probably  the  chief  cause  of 
the  dye  emanating  from  the  grooves  at  that  location.  The  blemish  can  be 
considered  an  undesirable  obstruction  that  slows  and  ponds  the  runoff.  Figure 


13(b)  Is  more  typical  of  the  results  of  this  research  A  closeup  photograph 
of  a  more  uniform  section  of  the  slab  shows,  on  the  right,  dye  diffusing  out 
of  the  grooves  because  they  are  filled  by  all  of  the  Inflow  upstream.  Motion 
pictures  of  the  experiments  were  also  taken,  and  are  in  the  posession  of  FAA 
personnel  at  the  Technical  Center  In  New  Jersey. 

The  most  significant  observation  from  these  dye  tests  was  that  all  the 
rain  falling  on  the  upstream  end  was  carried  in  the  grooves  in  such  a  way  that 
the  slab  surface,  although  wet,  had  zero  water  depth.  The  width  of  the  area 
contributing  lateral  Inflow  to  the  grooves  would  be  equal  to  the  groove 
spacing  for  spaclngs  of  5  In.  and  less.  At  what  larger  spacing  this  condition 
would  cease  to  occur  is  uncertain.  By  programming  this  condition  into  the 
computer  model.  It  will  be  possible  to  predict  the  location  where  the  grooves 
begin  to  overflow  and  surface  depths  again  start  to  build  up. 


4.  ANALYSIS  OF  RESULTS 


4.1  MODIFICATIONS  TO  THE  COMPUTER  MODEL  OF  PHASE  I 

The  FORTRAN  IV  computer  program,  discussed  In  Section  2.3,  was  modified 
to  reflect  the  experience  gained  from  the  experiments  described  In  the 
previous  chapter.  The  Manning's  n  value  used  In  Phase  I  for  the  pavement 
surface  was  lower  than  the  experiments  have  Indicated,  but  since  It  was  an 
Input  value  to  the  original  program,  no  program  modification  was  necessary. 
However,  the  program  was  modified  to  accept  as  Input  a  separate  Manning's  n 
value  (ng)  for  grooves.  Another  modification  to  the  program  was  based  on 
the  concept  that  all  of  the  rain  falling  at  the  upstream  end  of  the  grooved 
pavement  Is  carried  In  the  grooves  and  not  on  the  surface.  This  condition 
continues  until  the  grooves  fill  up  at  some  point  downstream.  At  that  point, 
the  pavement  surface  water  depths  begin  to  be  computed  as  though  flow  were 
occurring  over  a  plane  surface  with  constant  n.  The  program  assumes  that  all 
rain  falling  on  a  spacing  of  5  In.  or  less  contributes  lateral  Inflow  to  a 
groove.  The  Phase  I  computer  program  used  the  wetted  perimeter  of  a  groove  as 
an  arbitrary  determination  of  the  rainfall  contributing  to  flow  in  the  groove. 
While  the  wetted  perimeter  concept  remained  the  basis  of  the  model 
quantitatively,  the  original  computer  program  was  modified  to  use  the  actual 
geometry  of  grooved  pavement  because  there  was  no  flow  on  the  upstream 
pavement  surface.  The  nature  of  this  model  made  adjustment  to  actual 
conditions  easy  by  Incrementally  comparing  rain  volume  rates  to  groove  flow 
downstream  to  determine  the  point  where  the  grooves  overflow.  The  origin  of  L 
In  equation  (9)  Is  then  shifted  to  this  point  for  subsequent  computations  of 
water  depths  on  the  pavement  surface. 

The  flow  diagram,  listing  of  statements,  guide  to  Input,  and  typical 
output  of  the  modified  computer  program  are  presented  In  Appendices  C,  D,  E, 
and  F,  respectively.  The  program  was  executed  on  an  IBM  370/Model  3033 
system.  The  modified  model  results  are  shown  graphically  In  this  chapter. 

The  Manning's  n  values  shown  on  graphs  from  this  point  on  are  those  for  the 
pavement  surface  and  not  for  the  grooves.  The  groove  value,  ng  s  0.010,  was 
constant  for  all  of  the  computer  runs.  Reflex-percussive  grooves  were  not 
Included  In  these  runs  because  their  ng  value  Is  unknown,  and  because  they 
were  not  Included  In  the  experiments. 


4.2  EXPERIMENTAL  VERSUS  ANALYTICAL  WATER  DEPTHS 


Figures  14  through  18  display  comparisons  of  experimental  water  depths 
and  analytical  curves  of  water  depths  generated  by  the  modified  computer 
program.  Each  of  the  five  graphs  represents  a  specific  combination  of  one  of 
the  two  rainfall  Intensities  with  either  the  ungrooved  or  the  grooved  slab  at 
a  particular  spacing. 

Figures  14  and  15  for  the  ungrooved  pavement  show  three  computer  model 
curves,  one  for  each  of  the  three  pavement  surface  n  values  discussed  in 
Section  3.3.1.  The  best  comparison  of  data  and  model  appears,  surprisingly, 
to  be  for  n  ■  0.081,  which  Is  the  upper  limit  value  from  Figure  4  using  the 
average  texture  of  the  pavement  blemish.  It  would  be  difficult  to  conclude, 
however,  that  the  blemish  controls  the  hydraulic  resistance  of  the  entire 
slab.  Figures  16,  17,  and  18  clearly  show  the  reduction  of  water  depths  as  a 
result  of  grooving.  Ir  Figure  18  the  experimental  depths  are  so  small  that 
equation  (11)  might  no:  be  satisfied.  Hence,  the  Manning  equation  used  In  the 
kinematic  wave  approximation  In  the  computer  model  may  not  be  valid.  This  may 
account  for  the  poor  comparison  of  data  and  model. 

Two  types  of  errors  occurred  In  this  study:  errors  which  only  affect  the 
measurement  of  depths  and  those  which  affect  the  comparison  of  the 
experimental  results  with  the  computer  model  results.  The  possible  sources  of 
error  are  defined  here,  but  their  magnitudes  are  uncertain. 

Several  possibilities  exist  for  depth  measurement  errors.  The  flexible 
tubes  that  extend  to  the  sensor  openings  In  the  slab  were  sealed  to  the  slab 
with  a  commercial  windshield  sealer  compound.  The  holes  are  so  small  that 
uniformity  In  sealing  was  not  possible.  In  some  cases,  the  sealer  was  rough 
and  bulged  up  slightly  around  the  top  edges  of  the  tubing.  The  zero  water 
depth  datum  for  each  sensor  opening  Is  set  according  to  the  judgment  of  the 
researcher.  The  effect  of  capillary  action  on  this  setting  is  uncertain,  as 
Is  the  effect  of  using  a  mild  detergent  as  a  wetting  agent  to  control  this 
condition  by  reduction  of  surface  tension.  Perhaps  a  better  solution  in 
future  studies  would  be  to  countersink  the  sensor  openings  on  the  slab  prior 
to  placing  and  sealing  the  tubing.  Another  source  of  error  Is  the  possibility 
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gure  18.  Water  Depth  versus  Distance  for  2.5-inch  Rectangular-Groove  Spacing  at  2.45 
inches  per  hour 


of  air  in  the  tubing,  which  could  cause  a  false  transmission  of  water  depth 
pressure  to  a  transducer.  Since  the  tubing  beneath  the  slab  Is  not  visible, 

not  only  is  this  condition  difficult  to  control,  it  is  difficult  even  to 
check.  A  system  where  the  transducer  is  located  Immediately  beneath  the  slab 
would  probably  be  more  reliable.  In  addition,  the  calibration  of  each 
transducer  was  accomplished  very  well,  but  not  exactly  to  a  100*  correlation, 
to  which  Figure  12  will  attest.  Then,  too,  the  fluctuation  of  voltage 
readings  on  the  digital  multimeter  required  judgment  on  the  part  of  the 
researcher  to  record  a  single  reading.  Such  fluctuations  In  voltage  might  be 
augmented  according  to  whether  a  sensor  opening  was  subject  to  periodic  Impact 
from  pelting  raindrops  or  a  transducer  was  kept  absolutely  dry.  Also,  In 
repeated  tests  the  ability  to  create  exactly  the  same  rate  and  pattern  of 
rainfall  by  setting  the  simulator  equipment  Is  uncertain.  Furthermore, 
fluctuations  in  voltage  to  the  pump  supplying  the  rain  were  never  the  same 
from  one  run  to  another.  If  arbitrary  values  are  given  to  the  sources  of 
errors,  It  Is  estimated  that  any  depth  measurement  could  be  subject  to  an 
error  of  perhaps  ±0.04  In. 

In  comoaring  the  model  results  with  the  measured  water  depths, 
differences  between  what  the  computer  assumes  In  Its  solution  and  what 
actually  occurs  In  the  experimental  system  become  sources  of  error.  The 
computer  model  assumes  an  absolutely  perfect  plane  sloping  at  1.5*  with 
texture  depth.  Manning's  n,  and  rainfall  intensity  perfectly  uniform.  The 
slab  was  constructed  within  forms  set  accurately  to  1.5*  slope.  The  surface 
finishing  of  the  concrete  was  accomplished  as  carefully  as  cost  would  permit, 
but  some  slight  ponding  was  evident  during  tests.  The  blemish  Is  the  primary 
example  of  the  texture  not  being  uniform.  Also,  the  slab  may  have  settled 
slightly  since  construction,  perhaps  nonuniformly,  since  a  small  tension  crack 
is  visible  on  the  surface.  The  rainfall  simulation  equipment  does  not 
generate  absolutely  uniform  rainfall  rates:  there  are  spatial  variations  in 
rainfa.i  intensities.  In  the  development  of  the  equipment  these  variations 
were  minimized.  Still,  standard  deviations  of  the  spatial  average  intensities 
in  this  study  may  be  of  the  order  of  ±0.3  in. /hr.  The  spatial  variations 
would  affect  not  only  the  depths  down  the  slab,  but  also  the  lateral  average 
deoths  shown  In  the  Appendix  B  data.  Certainly  the  experimental  errors  in 
measurinq  depths  also  affect  the  comparison. 


It  does  not  seem  logical  to  absorb  all  the  errors  of  experimentation  into 
an  adjusted  value  of  Manning's  n,  but  numerous  other  studies  have  done  so. 

The  question  is  whether  an  n  value  can  be  accurately  predicted  for  a  given 
situation,  or  whether  It  must  be  measured. 

4.3  PREDICTION  OF  HATER  DEPTHS  FOR  RUNWAY  LANES 


The  modified  computer  program  was  used  to  predict  equilibria  water 
depths  for  rectangular  grooves  at  the  four  spaclngs  out  to  a  width  of  runway 
lane  of  100  ft.  To  justify  the  selection  of  a  Manning's  n  value  for  these 
computer  runs,  the  upper  limit  from  Figure  4  was  used  with  an  average  texture 
depth  of  0.026  In.  to  obtain  n  ■  0.056.  The  rainfall  rate  selected  was 
3  in. /hr,  the  next  highest  Integer  value  beyond  the  maximum  used  In  the 
experiments  and  near  the  maximum  achievable  by  the  rain  simulator.  The  output 
for  these  computer  runs  Is  shown  in  Appendix  F.  The  runs  are  summarized 
graphically  In  Figure  19. 

It  appears  on  the  graph  In  Figure  19,  and  can  be  shown  by  equations, 
that  there  Is  an  inverse  linear  relation  between  spacing  and  the  value  at 
which  each  curve  Intersects  the  distance  axis  (the  predicted  point  where 
grooves  fill  and  overflow).  If  the  distance  coordinate  were  multiplied  by  the 
spacing,  the  curves  would  collapse  Into  one.  This  has  not  been  done  since  it 
was  judged  that  the  relation  among  the  variables  Is  clearer  as  depicted.  It 
should  be  noted  that,  on  the  Appendix  F  output  sheets,  equillbriun  time  at  the 
100-ft  edge  Is  the  same  for  all  four  grooving  conditions.  This  assumption  is 
reasonable  since  the  equilibrium  flow  rate  off  the  edge  of  the  pavement  or  at 
any  other  given  distance  downstream  must  be  the  same  for  all  four  cases  at 
3.00  In. /hr.  The  principle  of  conservation  of  mass  would  be  true  for  any 
hydraulic  roughness  as  well,  although  hydraulic  roughness  was  not  varied  for 
the  four  cases.  The  following  calculation  will  Illustrate  that  total  flow  off 
the  edge  of  the  runway  must  be  equal  to  total  flow  onto  the  runway  at 
equilibrium  conditions: 


Flow  onto  pavement  ■  3  In. /hr  x  1  ft/12  in.  x  1  hr/3600  sec 

x  100  ft  x  15  ft 
■  0.104  cfs 


A  check  of  the  Appendix  F  output  will  confirm  that  this  figure  agrees  with  the 
flow  off  the  edge  of  the  runway  as  calculated  from  pavement  roughness  and 
grooving  considerations  by  using  equation  (10). 

From  Figure  19,  it  can  be  observed  that  the  predicted  maximun  reduction 
in  water  depths  due  to  grooving  Is  about  0.18  In.  for  rectangular  grooves 
spaced  at  1.25  In.  at  a  distance  of  about  39  ft.  The  same  point  is  shown  on 
Figure  20  as  a  100%  reduction  In  depth.  The  curves  of  Figure  20  also  collapse 
into  a  single  curve  in  a  manner  similar  to  those  In  Figure  19. 

Figure  21  represents  curves  of  prediction  of  water  depths  for  the 
reflex-percussive  grooves  of  Figure  1  at  the  same  surface  roughness  (0.056) 
and  rainfall  Intensity  (3.00  In. /hr)  as  for  the  rectangular  grooves  in  Figure 
19.  However,  the  groove  roughness  (ng)  is  uncertain  for  these  grooves.  The 
points  necessary  to  plot  these  curves  were  computed  manually,  since  no 
provision  was  made  for  this  shape  In  the  modified  computer  program.  It 
appears  that  at  the  same  ng  value  of  0.010,  these  grooves  are  much  more 
efficient  In  reducing  water  depths  than  are  rectangular  grooves  spaced  at  2.5 
in.,  and  almost  as  efficient  as  rectangular  grooves  spaced  at  1.25  in. 


Predicted  Percent  Reduction  in  Depth  versus  Distance  fcr  Various  Rectangular 
Groove  Spacings  at  3.00  inches  per  hour 


Figure  21.  Predic 


5 .  SUMMARY 


The  subject  of  this  study  has  been  the  prediction  of  runoff 
characteristics  from  a  laterally  grooved  runway  during  a  steady,  uniform 
rainfall.  The  conclusions  and  recommendations  will  address  this  situation. 
However,  the  authors  believe  that  one  of  the  principal  advantages  of  grooving 
Is  the  quick  draining  of  the  pavement  after  rain  has  stopped,  a  point  also 
made  by  Horne  [26]. 

5.1  CONCLUSIONS 


1.  The  kinematic  wave  approximation  to  the  shallow  water  equations, 
together  with  an  estimate  of  hydraulic  roughness  (Manning's  n),  is  a 
reasonable  model  for  predicting  runoff  characteristics  from  a  grooved 
or  ungrooved  runway. 

2.  The  lateral  grooving  of  runways  enhances  drainage  In  the  form  of 
decreased  water  depths.  The  following  table  gives  a  quantitative 
Indication  of  the  decreases  predicted  from  Figures  19  and  20  for  3.00 
Inches  of  rain  per  hour  and  n  ■  0.056. 


Groove 

Spacing 

(in-) 

100%  Reduction  in 

Water  Depth  through 

L  (ft) 

Maximum  Reduction 

In  Water  Depth 
(in.) 

%  Reduction 

In  Water  Depth 

at  L  *  100  ft 

1.25 

39 

0.18 

27 

2.50 

19 

0.11 

14 

5.00 

9 

0.07 

7 

3.  The  efficiency  of  a  groove  shape  in  reducing  water  depths  Is  directly 
dependent  on  Its  volumetric  capacity  and  inversely  dependent  on  Its 
hydraulic  roughness  coefficient  (Manning's  n). 
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4.  The  topography  of  a  runway  surface  subject  to  steady,  uniform  rain 
will  determine  the  magnitude  of  equilibrium  runoff  depths  that  will 
occur  there;  however,  it  will  not  determine  the  magnitude  of 
equilibrium  runoff  flow  rates,  which  must  obey  the  law  of 
conservation  of  mass. 

5.2  RECOMMENDATIONS 


A  logical  continuation  of  this  research  would  be  to  repeat  the 
experimental  analyses  for  an  asphaltic  concrete  surface  (probably  a  wearing 
course)  laid  over  the  existing  Portland  cement  concrete  slab.  The  effects  of 
an  adverse  wind,  and  possibly  a  variable  rainfall  rate,  on  the  pavement 
drainage  should  also  be  studied.  By  retaining  the  groove  shape  as 
rectangular,  a  commonly  occurring  field  condition  would  be  simulated.  Such  a 
study  might  be  conducted  In  12  to  15  months.  Looking  even  further  ahead,  the 
same  runoff  tests  should  be  conducted  for  asphaltic  concrete  surfaces  with 
reflex-percussive  grooves.  Depending  on  what  parameters  are  chosen,  a  study 
of  this  nature  might  be  completed  In  an  additional  12  months. 

Also,  additional  experimental  and/or  theoretical  work  should  be  performed 
as  follows: 

1.  With  reflex-percussive  grooves,  particularly  with  regard  to  their 
hydraulic  roughness; 

2.  To  refine  further  the  notion  of  hydraulic  roughness  versus 
macrotexture,  as  depicted  In  Figure  4,  by  carefully  controlled 
channel  experiments  with  typical  runway  pavement  surfaces; 

3.  On  the  effects  of  ponding  on  runways,  which  may  be  more  crucial  than 
runoff  depths  from  the  standpoint  of  hydroplaning  (the  technology  to 
determine  accurate  runway  topography  would  have  to  be  developed  to 
detect  the  slightest  rutting,  and  perhaps  a  two-dimensional  runoff 
model  could  be  developed); 

4.  On  how  the  flow  occurs  over  grooves  once  they  are  filled.  This 
research  assumed  a  planar  flow  over  the  filled  grooves  at  an  n  equal 
to  that  of  the  pavement  surface. 
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APPENDIX  A.  TECHNICAL  DATA  ON  PRESSURE  TRANSDUCERS 


MANUFACTURER:  Foxboro/ICT  Inc.,  San  Jose,  California 
DESCRIPTION:  Series  1813-003-1 

A  linear  integrated  circuit  transducer  which  senses  pressure 
by  means  of  a  Wheatstone  bridge  diffused  Into  a  silicon 
diaphragm.  The  output  signal  Is  linear  and  proportional  to  the 
applied  pressure. 

PRESSURE  RANGE:  0  to  10  in.  of  H20,  gage 
INPUT  SIGNAL:  1.5  mA  (DC) 

OUTPUT  SIGNAL:  15  mV  ±  5  mV 
ZERO  READING:  0  ±  5  mV 
RESPONSE  EQUATION:  Y  -  K  (AV) 


where 


Y  ■  change  In  pressure  head  (In.  of  H20) 
aV  3  change  In  voltage  (mV) 

K  ■  calibration  coefficient.  Individual  for  each 
transducer  (In. /mV) 


CALIBRATION  COEFFICIENTS  (K)  (also  see  Figure  6): 


TRANSDUCER 


0.056 

0.063 

0.063 

0.068 

0.074 

0.074 

0.083 

0.083 

0.090 

0.094 


7rSWw??i!f 
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APPENDIX  B.  EQUILIBRIUM  WATER  DEPTH  DATA 


Table  B-l.  Average  Water  Depths  (in.)  at  a  Rainfall  Intensity 
of  1.15  in. /hr  on  an  Ungrooved  Surface 


01  stance 
(ft) 

A 

Location 

BCD 

Lateral 

Average 

15.6 

(Range) 

0.068 

0.064-0.075 

-- 

0.068 

20.7 

0.082 

0.070-0.098 

— 

0.079 

0.071-0.095 

0.080 

24.8 

0.093 

0.085-0.104 

0.108 

0.097-0.112 

0.088 

0.080-0.099 

0.095 

29.5 

0.169 

0.142-0.212 

0.118 

0.097-0.128 

0.121 

0.115-0.142 

0.136 

*  Blanks 

Indicate  data  not  taken  or  unreliable  readings. 

Table  B-2.  Average  Water  Depths  (in.)  at  a  Rainfall  Intensity 
of  2.45  In. /hr  on  an  Ungrooved  Surface 

Distance 

Location 

Lateral 

(ft) 

A 

B 

C  D 

Average 

15.6 

(Range) 

— 

— 

— 

— 

20.7 

0.132 

0.112-0.147 

— 

0.118 

0.109-0.133 

0.125 

24.8 

0.142 

0.122-0.155 

— 

0.152 

0.120-0.164 

0.147 

29.5 

0.281 

0.223-0.302 

0.191 

0.180-0.229 

0.185 

0.174-0.208 

0.212 
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Table  B-3.  Average  Water  Depths  (in.)  at  a  Rainfall  Intensity  of 
1.15  in. /hr  on  a  Grooved  Surface  with  5-in.  Spacing 


Distance 

(ft) 


Location 


Lateral 

Average 


15.6 

(Range) 

•• 

0.000 

20.7 

0.000 

0.000 

0.000 

— 

0.000 

24.8 

0.000 

— 

0.000 

— 

0.000 

29.5 

0.052 

0.039 

__ 

0.046 

0.048-0.069  0.033-0.058 


Table  B-4.  Average  Water  Depths  (In.)  at  a  Rainfall  Intensity  of 
2.45  In. /hr.  on  a  Grooved  Surface  with  5-1n.  Spacing 


Distance 

(ft) 


Location 


Lateral 

Average 


E  15.6 

i  (Range) 

*»«» 

0.040 

0.031-0.045 

0.040 

j  20.7 

0.075 

0.060-0.081 

0.071 

0.053-0.097 

-- 

0.073 

j  24.8 

0.118 

0.103-0.128 

0.109 

0.092-0.125 

0.102 

0.095-0.113 

0.110 

29.5 

0.159 

0.142 

0.138 

0.146 

0.131-0.187 

0.131-0.158 

0.129-0.155 

Table  8-5.  Average  Water  Depths  (In.)  at  a  Rainfall  Intensity  of  2.45  in. /hr 
on  a  Grooved  Surface  with  2. 5- in.  Spacing 


Distance 


Location 


Lateral 


APPENDIX  C.  FLOW  CHART  FOR  MODIFIED  COMPUTER  PROGRAM 


Input  Data 


Convert  Oata  to  Proper 
Units 


Computer  aX 


Is  Pavement  Grooved? 


Compute  tne  Following  Values 

•  Distance  Along  Flow  Path 

•  Equilibrium  Time  at  End 
of  Segment 

•  Equilibrium  Depth  at  End 
of  Segment 

•  Equilibrium  Flow  for 
Total  Plane 


Is  Groove  Length  Exhausted? 


Compute  the  Following  Values 

•  Groove  Flow  Capacity 

•  Distance  to  Reach  Groove 
Flow  Capacity 

•  Equilibrium  Time  to  Reach 
Groove  Flow  Capacity 


Compute  Distance  Along  Flow  Path 


Is  Distance  Less  than  Olstance 
to  Reach  Groove  Capacity? 


Is  Path  Length  Exhausted? 


Compute  New  AX 


Compute  the  Following  Values 

•  Olstance  Along  Flow  Path 

•  Equilibrium  Time  at  End 
of  Segment 

•  Equilibrium  Depth  at  End 
of  Segment 

•  Equilibrium  Flow  for 
Total  Plane 


Compute  the  Following  Values 

•  Equilibria  Flow  per  Groove 

•  Equilibrium  Flow  for 
Total  Plane 

•  Ratio  of  Flow  per  Groove  to 
Groove  Flow  Capacity 


Set  Equilibrium  Flow  Depth 
for  Section  Equal  to  Zero 


|  Output 


-v»y. 


*  ■ vWW.nV vW'.WW' 
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APPENDIX  0. 


LISTING  OF  MODIFIED  COMPUTER 
PROGRAM  IN  FORTRAN  IV 


8-0  (0  V0A/R1--1***) 


Tits  P1PCRA1  IS  DESICRED  TO  PRIDIC?  I|TU  DEPTHS  OH  1  PL  IKE  SQsrACE. 
VARIABLES  XHCLODE  S01P1CR  TUTOR?:,  SLOPE,  RAt If ALL  IBTESSITT  AX 0 
BAIHXHOS  ROCGIUISS.  ALLOP  ABC  I  XS  BASE  TOR  ONGROOTED  PAVEREBT  OR 
■ECTARCOLAR  GROOVES.  TEE  OROOVSS  CAR  BE  Or  AST  SIZE,  SPACIRO  OR 
LERCTR. 

BEAL  f  EQ(200)  ,  TPL  ABE  (20 3 »  ,  ttRCB  (2)0)  ,QPL  (200)  .BIST  ( 
f  240)  ,QTOT  (200)  ,  PIO¥S<200)  ,GPATIC(200) 

READ,  J 

DO  101  K  •  1.J.1 
PRIRT  1*0, R 
1*0  PORRAT  ('I*,  12) 

PRIRT,  •  • 

PRIRT, • 

SAB  OOTPOT' 

PBXRT,  • 

SVERSITf ' 

PRIRT,  • 

SR  XB9TIT0TS' 

PRIRT,  • 


READ,  TO, BOnOR 1 , ROBOBO 
BEAD,  BP 
READ,  SLOPE 

BRAD,  ISRAPE,GDtPTH,GBIDTH»0 SPACE 
BEAD,  PLERC, GLtRC.R ICT  R 


IP  (ISHAPS  .EO»  0)  PRIRT.  •  ORGRCCVED  SOPTACE* 

IP  (TSBAPS  .BO-  I)  PRIRT,  '  BECTISCtlLA*  GROOVES' 
PRIBT.’  • 

CORPERSIOR  OP  1PPUT  CATA  TO  POOT-S’COXD  OBITS 

QBP  -  BP/B3200 

SPACE  »  eSPACE/12 

BBS  ■  0VIDTR/12 

8RD  •  0DEPTR/1  2 

DELTA!  ■  Olf RO/100 

ALPBA  -  ( 1 .  BB/IO03BP)  •  (SLCPE)  ••  (0 . *) 

MOOTS  CORPI  CO  RATIOS  CETtPRIRATIOR 

IP  (XSRAPS  .  EO.  I)  0°  TO  10 

DEPTH  CALCRLATIORS  ORGROCV JO  CARP 

DO  10  2  I  *1,100,  * 

TEO(X)  •  1/QRM  (DELTA!  •  X  •  CEP/ALPSA)  ••  (0.*0) 
TPLARE(X)  •  CEP'TIQCI) 

QPl(X)  «  AlPflA  •  ((TPLIREd)  )*>(1.»*T))*SPACE 
TXRCR(t)  •  fPLAR  P  (I)  *12 
OTOT(I)  •  OPLIII  •  (RIDTB/SPACX) 

DIRT  (7)  •  DELTU*  X 
102  COBIXBOE 
BC  TO  *0 


c 

c 


OtPTR  CALCULATIONS  GAOOTCO  CAST 


10  rtOM*  1.  »7/iOrjCHO*CiD*G|s*  (C(N*GPD/  (GRU*7*gAD)  )  ••  (0.  66?)  * 
aslope**<o.S) 
ao  TO  30 

25  Mill,  •  Invalid  ji«»  *aaahete*  • 

80  TO  777 

JO  IP  (1 5  SAPS  .RE.  0|  QIBNO.rl.ONG/  (Off'SPACA) 

I  •  1 

*0  IP  ( (DELTA!*!)  .  Lt.  6(C«n)  THEN  00 

TEO(I)  ■  1/QAP*  (DBlTAX»I*aiP/AlP«l)  ••(0.60) 

Q»l  (II  •4lP«A»(0*r»TlcfI))*»  ( 1 . 6667) *SPACI 

OIATZO(T)  ■  0PL(X)/PLO«Q 

TIRCN(I)  -0.0 

DZST(I)  •  DELTA!  •  X 

QtOT(t)  *QPl(I)  *RICTP/SPACE 

I  •  1*1 

80  TO  AO 

ELSE  00  1 

TEQCP*  1/QA(*(GR INC *QAf/ ALPHA)**  (0.60) 

L*I 

DSTSP  •  DELTA!  *(L)  -  6RIR0 
EHO  IP 

00  *5  I  •  1,100, 1 

TB-Jd)  *  </Q*P*  ( (OPLTUd-CEESD)  •QRP/ALIRA)  ••(0.60) 
TPLAAE(I)  •  C*l«T*0(I) 

TIRCH  (I)  -  TPLARE(t)  •  12 

QPLCI  •  ALPHA  •  (  (T7LANR  (T) )  ••  ( 1 . 667 ) )  *3P  ACE  ♦  PLOVC 
OTOT(I)  •  OPI.(I)  •HXCT6/SPACP 
01  ST  (I)  •DELTA!*! 

CRATIO(I)  -A. 5 
*5  CONTINUE 

60  IP  (I SHAPE  .(0.0)  THEN  00 
DC  96  11-1,200,  1 

GfeATIO  |S)  -6.5 
99  CONTINUE 

teqi-.a-o.o 

ELSE  DO 
■  •  0 
1»0  IP 

OEPTH  CALCULATIONS  TOR  CASE  HHESE  CPOOTPS  00  ROT  PUN 
ENTIRE  PAT ZSENT  LIRCIB 

IP  (PLERC  .GT.  CLING  )  TRET  CO 

TPLARS  (  109)  -  (QPL  (100)  /  (ALPHA  ‘SPACE)  )  ••  (0.6) 

DX  *  (PLE.TS-CLENGJ/10.0 
TIRCH  (107)  •  IPLARE  ( >33)  •  12 
00  50  I  •  101,110,1 

TBOm-  1/0RP*(»0X *(1-110)  *0(ITAX*100)*QSP/ALPHA)**  (0 
TPLARB(I)  •  Q»r*TBC(I) 

QPL(I1  •  ALPHA  •  (  (TPLARC(l)l**  (1.667)  )  *SPAC» 

TIRCH  (I)  •  TPLAHIU)  *12 
QTOT(t)  *0?L|I)  •  (NIDT1I/SPAC*) 

6PATI0  (I  |  *  i*  .5 

OUT  (!)  •  0X*(I>103)  *  DELTA!*  100 
SO  CONTINUE 

TIN2C0  *  1/CAP*(PlfTG*'JPP/ALP*A|  **(3.63) 

PLOHE’J  •  0707(110) 


AWV.'.VX 


I (NO  •  110 
USE  DC 

TIREEQ  •  VQAP*(PLENfi*QAP/ALPHA)»*(O.SO) 
PLONSQ  •  OTOT(IOO) 

I  END  •  100 

IP  (ISHAPS  .IQ.  0)  6REND*0 .0 
EDO  IP 


O0TP9T  RC1TI  NE 


HUT, 

PRIRT, 

PRIRT. 

PRINT. 

PRIRT. 

PRIRT. 

PRIRT, 

PRIRT, 

PRIRT. 

PRIRT. 

PRIRT. 


TEXTURE  0IPTR  (INCHES):  '  ,TD 

fl  ANN  IRC  S  "N"  P  ACTOR  (SOP  P  ICE)  :*,POOCHP 
HANNINGS  *«•  P ACTOR (GROOVES) : ' .R09GHG 
RAIRPALL  1RTESSITT  (IN/RR)  t  '.  HP 
RORRAT  SLOPE  (PT/PT)  :  SLOPE 


6R00TI  SPACING  (IN):  *.  GSPACE 
OROOTE  DEPTH  (IN):  •,  CDEPTR 
GROOVE  RICIH  (IN):  •.  GVIOTR 


SORPACZ  LENGTH  (PT) 


IP  (ISHAPf  .  EO.  2)  Gtl*C*0.0 


PRIRT. 

PRIRT. 

PRIRT. 


OROOTE  LENGTH  (PT)  1  ■  ,  CLING 
SORPACE  NIOTH  (PT)  l',  RIDTH 


PRIRT.  ■ 
.PRIST.  • 
PRIRT,  • 
PRI*T.  • 
PRIRT.  • 
PRIRT.  • 


TOTAL  PLOT  AT  IR9  OP  PLANE(CPS)  :  *,  PLOREO 
TOTAL  SQOILIRPIOH  PLOR  TIRE  (SEC):  TIHEEQ 

GROOVE  EQOXLIRRIOS  LENGTH  (PT) :  '.GRERD 


PRIRT.  •  SEGHIJfT  DISTANCE  (PT)  PUN  D 

I  PLOR  (CPS /SPACING)  TOTAL  PLOV  (CPS)  GROOVE  PLOT 

PRIRT,  •  • 

DO  205  1-1.9. 1 

PRINT,  I.  OIST(I),  PINCH  (I),  QPl(I),  QTOT(I),  GRATIO(I) 
205  CCHTIROR 

DO  20*  1-10, IEH0.5 

PRINT, I.  DTST(I).  VXRCK(t).  «H(I),  OTOT  (I)  ,  CRATIO(I) 
20*  CORTINDE 
GO  TO  itt 
77T  R  -  0 
101  CORTINDE 
STOP 
IRO 


PIL.N  DEPTH (IN) 
GROOVE  PLOV  RATIO' 


••VANNING**  DNREPESENCRD  STATSSENT 


25  OSED  IN  LIRE 


57  FOLLOWS  A  TRANSFER 


STATE SERTS  RXICOTEO-  «0«I 


CORE  DS AGE 


OBJECT  CODE*  «720  SITES. ARRAP  ARPA-  ««00  SITES, TOTAL  AREA  AVAILABLE-  1«0*««  BITES 


DI AGNCSTICS 


RONBRR  Of  ERRORS- 


1,  RUHBER  OP  EXTENSIONS* 


CORPILE  TIPS- 


0.07  SRC.EXSCnTIOR  TINE*  0.10  SRC, 


RATPIV  -  HAR  1900  V2L0 


APPENDIX  E.  GUIDE  TO  COMPUTER  PROGRAM  INPUT 


This  guide  Is  a  modification  of  the  one  which  appeared  in  Reference  [1]. 

All  input  to  this  program  is  free  format.  If  more  than  one  variable  appears  on 
a  card,  they  need  only  be  separated  by  commas. 

An  input  data  deck  consists  of  one  J-card,  and  sets  of  A-E  cards.  Each 
separate  groove  pattern  requires  five  cards  (A,  B,  C,  D,  E-- in  that  order). 

J  Card:  Indicates  number  of  grooving  patterns  In  each  run 

Variable  Value  Description 

J  Integer  Number  of  grooving  patterns  in  the  run 

A-E  Cards:  A  complete  set  of  A-E  cards  must  be  input  for  each  desired  groove 

pattern.  Each  E  card  is  followed  by  an  A  card,  except  at  the  end  of 
the  grooving  pattern.  The  final  E-card  Is  followed  by  a  /*  card 
signaling  the  end  of  the  program. 


A  Card:  Provides  surface  roughness  Information 


Variable 

Value 

Description 

TD 

Real 

Depth  of  macrotexture  (In.) 

ROUGHP 

Real 

Manning's  n  value  associated  with 
texture  depth  on  the  pavement  surface 

ROUGHG 

Real 

Manning's  n  value  for  surface  of  groove 

B  Card:  Provides  rainfall  Intensity  Information 

Variable  Value  Description 

RF  Real  Rainfall  Intensity  (in. /hr) 

C  Card:  Provides  slope  information 
Variable  Value 

SLOPE  Real 


0  Card:  Provides 

groove  geometry  and  spacing 

parameters  (see  Figure  3) 

Variable 

Value 

Description 

I  SHAPE 

1 

Indicates  rectangular  (saw-cut) 
grooves 

0 

Indicates  ungrooved  surface 

GOEPTH 

Real 

Groove  depth  (in.).  ’  Use  0.00 
for  uhgrooved  surface. 

GWIDTH 

Real 

Groove  top  width  (in.).  Use 
0.00  for  ungrooved  surface. 

GSPACE 

Real 

Center-to-center  spacing 
between  adjacent  grooves  (in.). 
Use  WIDTH  x  12  for  ungrooved 
surface. 

E  Card:  Provides 

pavement  size  parameters 

Variable 

Value 

Description 

PLENG 

Real 

Total  distance  downslope  to  be 
studied  (ft) 

GLENG 

Real 

Total  distance  grooves  extend 
downslope  (ft).  Cannot  exceed 
PLENG. 

WIDTH 

Real 

Total  width  across  slope  to  be 
studied  (ft) 
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APPENDIX  F.  TYPICAL  COMPUTER  PROGRAM  OUTPUT 


RAA  BONW  AT  DRAIEIGf  PRCGPA1  OUTPUT 
PE  R.NSTLVANI  4  STAI?  ONTVERSITI 
PENMSTLVA.'II  A  TEA  NSPOPT  ATICN  INSTITUTE 
RO*!i3FF  1982 

ORJ5CO*  .  SURFACE 

TEXTURE  CFPTH  (INCUTS)  : 

■AIRINGS  TACTOR  (SOR  PACE)  : 

BARRINGS  "N"  IACT09  (CROOVES)  : 

RAINFALL  IETENSITT  (IN/HS)  : 

rorsa*  slope  (rT/rT)  : 

GROOVE  SPACING  (IN):  183.000000 

GROOVE  DEPTtl  (IN)  :  0.0 

GROOVE  RIOTH  (IN):  0.0 

SOFFiCE  LENGTH  (FT)  :  100.000300 

GROOVE  LIRCTR  (FT)  :  100.000000 

SURFACE  RIDTe  (RT)  :  IS.OOCOOOO 

TOTAL  RIOS  AT  END  OF  PLANE (CFS)  :  0.10403BJ 

TOTAL  EQOILTBRIUfl  now  TINE  (SIC):  359.  350930 

GROOVE  EQUIL1PR IUN  LENGTH  (PT):  0.0 


SIG.1ERT 

DISTANCE  (FT) 

WATER  OEPTH(IR) 

PLOW  (CPS/SPACING) 

TOTAL  PLOW  (CPS) 

GROOVE  PLOT  RAT 

1 

1.0000003 

0.0188946 

0.0010394 

0.0010394 

4.5099930 

2 

2.0000003 

0.0286  388 

0.0020791 

0.0020791 

4. 5000000 

3 

3.0000000 

0.0385267 

0.0031190 

0.0031190 

4.5000000 

A 

4.0000307 

0.0434983 

0.0041589 

0.0041589 

4. 5007900 

5 

5.00  C0000 

0.0496271 

0.0051988 

0.0051998 

4.5009000 

6 

6.0000031 

0.0553640 

0.0062388 

0.0062388 

4.5000030 

7 

7.0000000 

0.06  97290 

0.0072788 

0.0072788 

4.5000000 

a 

8.0000000 

0. 3657947 

0.0083189 

0.0083199 

4. 5009000 

9 

9.0000000 

0.0706127 

0.0093 590 

0.0093590 

4.5090090 

10 

10.0009000 

0.0752206 

0.0103991 

0.0103991 

4.5000090 

15 

15.00  00000 

0. 0959383 

0.0155998 

0.0155998 

4.5000099 

20 

20.0000090 

0.  1140132 

0.0209910 

0.0208010 

4.5009000 

25 

25.00 00100 

0.  113347  1 

0.0260024 

0.0260024 

4.5000000 

30 

30.0000009 

0.1454  152 

0.0312040 

0.0312040 

4.5000000 

35 

35.0001009 

0.  1595364 

0.0164058 

0.0164058 

4. 5000000 

«0 

40.00  00000 

0.1728118 

0.0416078 

0.0416076 

4.5000000 

♦  5 

45.9000033 

9.  1854662 

0.0468099 

0.0468099 

4.509.7000 

50 

50.0000001 

0.1975692 

0.0520120 

0.0520129 

4.5000030 

55 

55.0000009 

0. 2091967 

0.0572143 

0.0572143 

4.5000000 

60 

60.0003001 

0. 2204382 

0.0624166 

0.0624166 

4. 5000000 

65 

65.00  00000 

0.2312517 

0.0676191 

0.0676191 

4. 5000000 

70 

70.3003301 

0. 7417665 

0.0728217 

0.0728217 

4. 5000090 

75 

75.0000000 

0. 2519846 

0.0780243 

0.0780241 

4.5000000 

80 

80.0003300 

0. 26  >9335 

0.0812270 

0.0612270 

4. 5000070 

85 

85.0000009 

0.2716  36  7 

0.0884297 

0.0884297 

4.5090000 

90 

90.9033391 

0.  28 1 1 14  2 

0.0936376 

0.0936326 

4.  5003  00  9 

95 

95.0000000 

0. 24  13832 

0.0988  355 

0.0989355 

4.5000000 

100 

100.313103 

3.  244*589 

0.1040384 

0.  1040394 

4.5000000 

0. 0260000 
0. 0560000 
0.0103000 
3.0000000 
0.01S0D30 


•  w—  «r~~  %r  _ 


V>  *JI  "jr^w* -v  -T  ^7^7  w y ~»ir 


1ICTIRG0L1P  GPCOVZS 


TRtTJSE  DERI H  f I SCH  ZS )  : 
mviirgs  «*"  pactor  (snaPACO  : 
SHRliCS  FACTOR  (GROOVES)  : 

tinrm  irtzssiti  rm/HR)  : 
ror • ei  slope  (Pr/rt)  : 


GBCCVE  SP»Ct»G  (IS) : 
GROOVE  DEPTH  <I«)  : 
GRCOVJ  RIOTS  C*): 

SOP  f ICR  LERC18  ( fT) 
CPCCVR  LESGTR  (fT)  : 
30BPICE  SICTP  (PI)  : 


0. 0260  301 
0.3560100 

o. oi aaooo 

3.0000000 

0.0150000 

5.0000(130 
0. 2*00000 
0.  25300  00 

100. 030  '00 
100. 000000 
15.0000300 


lOTil  PLOD  IT  R  »D  OP  PURR  (CPS)  : 
TOTIL  RQOIUPPIOS  PLOP  PISE  (SEC): 
GROOVE  ZQOXUtRtOS  LEROTH  (?T)  : 


0.  1999941 
359.  3509)0 
9.9*30513 


SEGRERT 

1 

2 

3 

« 

5 

t 

7 

• 

9 

13 

15 

20 

25 

30 

35 

*0 

*5 

50 

55 

to 

45 

30 

75 

90 

95 

90 

95 

100 


DI5TASC  5  (PT) 

1.0000000 
2.0000007 
3.0000000 
9.0003000 
5.0000000 
9.0003000 
7.0000300 
9.0000000 
9.0000003 
10.0000 300 
15.00  00003 
20.0000000 
25.0000007 
30.0000003 
35.0000003 
90.0000030 
95.0000000 
50.0000003 
55.0000003 
60.0000000 
65.00 00301 
70.0000003 
75.0000003 
90.000003) 
95.0000703 
90.00  30001 
95.00C0007 
100.000000 


WATER  OR  PT*  (IP) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0032905 

0.09  99917 

0  .  0754593 
0  .  096  1  41} 
0.1191992 
0.  1335'24 
0.1955692 
0.  159651  1 
0.  1729  490 
0. 1955971 
0.  1976°»7 
0 . 2093 175 
0.  2235245 
0.  2313*97 
0.  24  19761 
0.2523912 
0.  2620  375 
0. 2717392 

0 . 29 12 1 33 


PLOW  (CPS/SMC IRC.) 

0.00002*9 
0.0000579 
0.0000966 
0.0001 157 
0.0091446 
0.0001736 
0.0002025 
0.0002314 
0.0002604 
0.0302993 
0.0004339 
0.0005792 
0.0007227 
0.0009672 
0.0010H6 
0.0011561 
0.0013006 
0.00 19951 
0.0015896 
0 • 00 17341 
0.0018796 
0.0020231 
0.0021677 
0.0023122 
0.3024667 
0.0026012 
0.0027457 
0.0028903 


TOTAL  PLOW  (CPS)  GROOf!  PtO*  RATIO 


0.03  1031* 

0.0020929 

0.0031294 

0.0041659 

0.0052073 

0.0062498 

0.0072903 

0.0083318 

0.0093734 

0.0104165 

0.0156154 

0.0239156 

0.0260164 

0.0312176 

0.0364190 

0.0416206 

0.0469224 

0.0 520241 

0 , 0572264 

0.0624286 

0.0676309 

0.07 28331 

0.0780357 

0.0832383 

0.0884498 

0.0936435 

0.0999463 

0.  1040491 


0. 1005100 
0.  2010232 
0.  3015374 
0.  4320520 
0 . 502S667 
0.6030820 
0.7035986 
0. 804i 148 
0.9006317 
9.5000030 
4.5000039 
9.5000000 
9. 5000099 
9.5000000 
9.5090979 
9.5000909 
9. 5000030 
9.5000000 
9.5000 000 
9.5000303 
9. 500000 0 
9. 5000000 
9.5000009 
9. 5000033 
9.5000070 
9. 5000090 
9.5000009 
9.5090000 
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KCTIHGOLl?  CHOOSES 


TEITOHE  DEPTH  (I  SCHES)  :  0.0260039 

bissisgs  pictoh  (sonnc?) :  o.osaoooo 

himisgs  «h”  r*CTot  (citooves) :  0.0190000 

iltSf ILL  IHTIHSITT  (IH/HS) :  3.0003900 

■mill  SLOPE  (PT/fT)  :  0.0150000 


CHOOSE  SPICIIG  (IS):  2.5000000 
GPCOSf  DEPTH  (IS):  0.  2500000 
CICCPI  SECTS  (IS):  0.2500000 


SOPflCE  LIIGTH  (PT)  :  100.000000 
CHOOSE  LESGTH  (PT)  :  100.000000 
SU5  I  ICE  SIOT S  (FT)  :  15.0000000 


TOTIL  FLOS  IT  ISO  Of  PLINE(CPS)  :  0.10*0601 
TOT IL  EQOtltPHIOS  PLOV  TIHF  (SEC):  359.  3508)3 
CPCOSE  EQ0UI96IU8  LESCTH  (PT)  :  19.89*1090 


SCCSEHT  DISTISCE(PT)  WITH  DEPTH(IH)  PLOV  (CPS/Sf ICISG)  TOTIL  PLOVtCPS)  CHOOSE  PLOV  IITIO 


1 

1.0000300 

0.0 

2 

2.0000000 

0.0 

3 

3.0000003 

0.0 

• 

8.0000000 

0.3 

5 

5.0000000 

0.0 

I 

8.0000000 

0.0 

7 

7.0000000 

0.0 

8 

8.0000000 

0.0 

5 

9.0000000 

0.0 

10 

10.0000000 

0.0 

15 

15.0000000 

0.0 

20 

20.0300303 

0.0099121 

25 

25.0000000 

0.0502551 

30 

30.0000000 

0.0758975 

35 

35.0000000 

0.0963*91 

•  0 

80.0000303 

0. 1183750 

•  5 

85.0000000 

0. 1306  781 

50 

50.0000000 

0.  1857230 

55 

55.0300307 

0.  1577958 

80 

80.3000001 

0.  17J0881 

85 

85. OOCOOOO 

0.  1857280 

70 

70.0000003 

0.  1918202 

75 

75.0000000 

0. 209*  38  3 

80 

80.0003000 

0.  2236*16 

85 

85.0000000 

0.7318778 

90 

90.0000000 

0.  28  1*858 

95 

95.0000000 

0.2571987 

100 

130.300700 

0. 2*2i*’» 

0.00301*5 

0.0010*1* 

0.0502550 

0.0009289 

0.0020829 

0. 1005116 

0.000083* 

0.003124* 

0.1507697 

0.0009579 

0.09*1659 

0  .  2010  269 

0.0000723 

0.0052073 

0.7512830 

0.0000868 

0.0062*98 

0. 3015*10 

0.0001013 

0.9072903 

0.  3517992 

0.0001157 

0.0083318 

0. *020573 

0.0001302 

0.00937  39 

0. *523160 

0.30019*7 

0.010*1*9 

0. 59257*3 

0.0002170 

0.015622* 

0.7538671 

0.0002893 

0.0208330 

*.5009000 

0.0003616 

0.0260319 

8.5000000 

0.0099338 

0.0312322 

8.5003097 

0.0095060 

0.0164370 

8.5000000 

0.9305783 

0.0*163*2 

9.5-700030 

0.0006505 

0.0*68356 

9.5000000 

0.9007227 

0.0520372 

8.500000-1 

0.9007950 

0.0572390 

8.5000000 

0.9008612 

0.0626*10 

8. 5090900 

0.0009  395 

0.0676*30 

• . 5000009 

0.0310117 

0.0728*52 

8.5007000 

0.90108*0 

0.0780*75 

8.5000000 

0.0011562 

0.0832*98 

9.  5000033 

0.0012285 

0.088*523 

8.5003000 

0.0013038 

0.09  365*7 

8. 5000303 

0.90137  33 

0.0988575 

8. 5000000 

0.001**53 

0. 10*0601 

8.5000000 

6 


*i: :**r.ct»E  uf.i’otf  j 


TOT XI.  FLOH(CrS) 

0.0010*1* 

0.0020829 

0 . 00  31299 

0.00*1659 

0. 0052073 

0.0062*88 

0. 0072903 

0.0083310 

0.009373* 

0.010*1*9 

0.015622* 

0.020*300 

0.0260376 

0.0312*53 

0. 036*529 

0.0*16661 

0.0*69651 

0.052065* 

0.0572662 

0.062*673 

0.06766*7 

0.0728703 

0.0780722 

0.08327*1 

0.088*762 

0.093678* 

0.0989807 

0.10*3830 


GfOOtS  FLO*  FXTIC 

0.0251275 
0.  0502550 
0.07538*3 
0. 1005129 
0. 1256*15 
0. 150773* 

0.  1758996 
0.  2013296 
0.2261578 
0.2512970 
0.  776*33* 
0.5025812 
0.62822*9 
0. 75387*2 
0.8795261 
*.5000090 
*.5009000 
*.5000000 
*.5009090 
*. 5000099 
• . 500''093 
*. 5030903 
*.5000000 
*.5000003 
*.5009090 
*.5009000 
*.5000090 
*,  50-90393 


